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25 1. Introduction 



26 Positron Emission Tomography (PET) is a functional imaging technique 

27 of growing importance in medical diagnostics. Its powers lie in the ability 

28 to reveal biologically significant processes that can be used, for example, in 

29 cancer screening and in studying neurodegenerative diseases. Conventional 

30 PET detectors employ scintillating inorganic crystals [l| as the gamma ray 

31 detection media. While crystal-based PET systems perform adequately for 

32 many applications there is motivation for seeking improvements of resolu- 

33 tions in energy, position, and time response to improve image quality and 

34 increasing overall sensitivity. Liquid xenon (LXe) is another gamma ray de- 

35 tector technology applicable to high resolution PET which may result in 

36 improved performance and reduced noise in images due to superior energy 

37 resolution, true 3-dimensional position reconstruction^and the capability for 

38 determining the Compton scattering sequence 0, 0, Isj. Energy resolution 

39 of 7% (FWHM) has been reported in small LXe detector tests by combin- 

40 ing scintillation light and ionization charge measurements j6|. Measuring 

41 charge in a drift chamber has been shown to provide 3-D sub-millimeter spa- 

42 tial resolution 0, @] because electron diffusion is very small [9] . In addition, 

43 sub-ns timing resolution has been achieved by measuring the scintillation 

44 light fld\. Liquid xenon is also inexpensive compared to crystal detectors 

45 commonly used for PET. Liquid xenon PET systems have the potential to 

46 reduce detector contributions to PET to the level of intrinsic limitations 

47 due to positron range and non-colinearity of the emitted photons. 

48 This paper deals with the energy resolution obtained from light and 

49 charge signals observed in a small LXe prototype detector as well as an 

50 investigation of the components influencing it and the sources of uncertainty 

51 which may inform the design of future detectors for PET. 

52 2. Micro-PET Detector Design 

53 We have developed a concept for a micro-PET detector shown in Fig. [1] 

54 that takes advantage of all the high resolution capabilities of LXe gamma ray 

55 detectors. Scintillation light is measured by arrays of large area avalanche 

56 photodiodes (LAAPD), which have been found to work well in LXe 

57 Charge measurement is achieved by using a time projection chamber (TPC), 

58 an approach successfully demonstrated in [Tj . Photons entering the LXe pro- 

59 duce prompt scintillation light and ionization which drifts under an electric 

60 field applied between the cathode and the anode of the TPC. The anode 



module (not shown in Fig. [T]) consists of a shielding grid followed by an 
array of wires preceding the anode which is segmented into strips perpen- 
dicular to the wires. The electron signal induced on the wires and collected 
by the strips provides a two dimensional (x-y) position measurement of 
the charge. The third coordinate (z) is obtained by measuring the elec- 
tron drift time i.e. the difference between the time of the light flash and 
the electron arrival time on the anode. Since every interaction is precisely 
recorded, Compton scattering can be reconstructed giving information on 
the direction of each incoming photon providing the possibility to suppress 
accidental coincidences and scattering prior to reaching the detector. 




Figure 1: The LXe PET ring concept. Scintillation light and charge are measured in 
each of the 12 modules consisting of a LXe time projection chamber viewed by avalanche 
photodiodes. 

The expectations for performance under operating conditions for PET 
include sub-millimeter 3-D position resolution from charge, timing resolu- 
tion of < 1 ns from scintillation light, energy resolution < 10% (FWHM) 
combining light and charge signals, and the ability to reconstruct Compton 
scattering. Spatial location of events obtained from the prompt distributed 
light signals will be used to reduce the ambiguities of associating the scintil- 
lation light and charge at high levels of activity. A simulation of the imaging 



performance of this system will be presented in a future publication 
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3. Small Chamber Prototype 

3.1. Test Setup 

As an initial step in studying LXe detectors for PET, we constructed 
a small test chamber (27 cm^) for simultaneous measurements of light and 
charge. The test chamber is shown schematically in Fig. [2l 
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Figure 2: Schematic views of the small test chamber. The side view illustrates the drift 
direction between the cathode and anode, viewed by two LAAPDs immersed in the LXe. 
The top view shows the segmentation of the anode. 



Scintillation light was detected by two 1.6 cm diameter windowless 
LAAPDs (Advanced Photonics Inc. [l6|). The LAAPDs were located at 
the center of the drift region as shown in Fig. [2], 1.5 cm above the grid wires 
and 1.5 cm below the cathode. Charge was collected on a central 1 cm di- 
ameter electrode (Al) or on an outer electrode (A2). An electric drift field 
was applied between the cathode and a shielding grid separated by 3 cm. 
The electric field was formed by a field cage consisting of 9 wires with a 
spacing of 3 mm strung along the four walls of the chamber. The voltage 
was distibuted by lOOMfi resistors. The APDs were outside the field cage 
and the distance between the field cage wires and the APDs was 2mm. The 
shielding grid consisted of 0.1 mm dia. wires spaced by 3 mm located 3 mm 
from the anode charge collection plane. The electric field between the grid 
wires and the anode was set higher than the drift field ensuring that all 



the electrons pass through the grid 13|]. In order to study the influence of 



98 the drift field on quantities like charge and light production and the energy 

99 resolution several settings were used. With the grid at ground potential 

100 measurements were made with the negative cathode voltage set to 1, 3, 6, 

101 and 8 kV, with the respective anode voltages set to 300, 600, 1200 and 

102 1200 V. Photons of 511 keV emitted after annihilation of positrons from a 

103 ^^Na source with an activity of 9.61 • 10^ Bq and situated in a collimator with 

104 an opening angle of 2° positioned 30 cm away from the cathode entered the 

105 test chamber (along the z axis) through the cathode plane. The trigger was 

106 generated by selecting signals in coincidence of both APDs and an external 

107 Nal detector placed at a distance of 50 cm from the source observing the 

108 full energy of the other 511 keV photon from the positron annihilation. The 

109 probability to detect more than one event in the chamber at the same time 
no was less than 3%. The detector was operated at 15 psia and at temperatures 

111 between 168 and 169 K. Before inserting the liquid in the vessel holding the 

112 detector a bake-out in vaccuum at 7.6-lO^^T was performed for 6 days at 

113 60°C to clean the components. The purification of the xenon was done in 

114 the gas phase using two stages both with equipment from From NuPure 

115 Corporation pjs]: first, the heated getter (NuPure Omni 600) was used to 

116 remove H2O, O2, CO, H2, and N2 to sub-ppb levels followed by a room 

117 temperature getter (Eliminator 600 eg) to remove H2O, O2, CO, H2, and 

118 hydrocarbons to < 0.5 ppb. The lifetime of drifting ionization electrons 

119 was used to indicate successful operation of the purification as discussed 

120 below. 

121 3.2. Readout Electronics 

122 The two anodes segments and the grid wires (ganged together) were 

123 connected to charge- sensitive amplifiers followed by a 1 ps time constant 

124 RC-CR shaper. The amplifier was calibrated using a narrow pulse input 

125 charge with a precision of 5%. The amplifier outputs were fanned out into 

126 three branches: 

127 1. A constant fraction discriminator followed by a time-to-digital con- 

128 verter (TDC) CAEN model V1190B; 

129 2. A charge sensitive analog-to-digital converter: 12 bit QDC CAEN 

130 model V792 with gate adjusted to the drift time and pulse shape; and 



131 3. A 20 MHz sampling waveform digitizer VF48 [15 



To get absolute charge values, the digitized waveform measured with the 
VF48 was used for the analysis presented in sections H] and O Because the 



134 QDC had a better signal to noise ratio it was used to determine the energy 

135 resolution in section [H] which did not require absolute charge calibration. 

136 The other reason the QDC was not used for absolute values was due to the 

137 very short pulse used for calibration. The longer chamber signal would not 

138 have been fully integrated within the window set. 

139 The observed range of noise of the amplifier was 700-1100 electrons due 

140 to varying external sources of induced noise. To reach the optimal position 

141 resolution, a signal to noise ratio larger than 5 was desirable requiring the 

142 electronics noise to be kept below 1000 electrons equivalent noise charge 

143 (ENC). A typical signal was expected to be at least 10 000 e-, as long as 

144 the electron attachment during the drift was small. 

145 The LAAPD voltages were set so that their gains were 500 and each was 

146 connected to a current-sensitive preamplifier with a pulse width of 50 ns and 

147 10^ electrons ENC. The amplifier signal was split into 3 branches: 

148 1. Discriminator and TDC; 

149 2. QDC CAEN model V792 with a gate of 100 ns; and 

150 3. 1 GHz waveform digitizer CAEN model V1729. 

151 Solid angle calculations showed that 12% of the scintillation photons 

152 reached the LAAPDs when the gamma interaction took place in the center 

153 of the chamber. 

154 4. Charge Collection 

155 The grid wires shielded the anodes from the current induced by the 

156 drifting electrons. Once the ionization electrons passed the grid wires the 

157 signal on the anode started to build up with a pulse shape that was largely 

158 independent of the z position of the primary interaction although electro- 

159 static calculations showed that the current pulse shape depended on the x-y 

160 distance of the electron cloud from the individual grid wires. Furthermore, 

161 depending on the drift velocity (typically 0.15 to 0.21 cm/jis), the current 

162 induced on the anodes lasted 1.5 to 2 ps. 

163 The waveforms measured on the grid and on the two anodes provided 

164 information about the location of charge creation and Compton scattering if 

165 multiple charge pulses were observed. The time of charge arrival relative to 

166 the light signal gave information about the position along the drift direction 

167 of the electrons. 

168 Figure [3] shows examples of four charge waveform events recorded by 

169 the 20 MHz waveform digitizers chosen to illustrate several types of events. 




Figure 3: Example of waveforms with a 2 kV/cm drift field. Central anode Al (dashed 
line), peripheral anode A2 (solid fine), grid (dotted line). 



170 In the upper two plots the charge was created roughly at the same z posi- 

171 tion; in the left plot the signal of the central anode (Al) integrated to zero 

172 and the anode A2 collected the charge, whereas in the right hand plot the 

173 interaction deposited the full charge on the central anode Al. The lower 

174 left panel shows an event originating close to the cathode plane, resulting 

175 in a measured drift time of 15 jis. The charge was collected by the central 

176 anode. The bipolar shape of the grid signal is clearly visible. The waveform 

177 measured on the grid depends on the z position of the interaction and is 

178 also influenced by electrons collected by the grid. The lower right panel 

179 shows two photon interactions, presumably one Compton scattering and 
iBo one photo-electric interaction. One interaction took place above Al and 
iBi one above A2. The interactions were also separated in the drift direction so 
1B2 that a two peak structure is visible in the grid waveform. Simulations of the 
183 setup showed that only a small fraction (less than 5%) of events fully con- 



1B4 tained on Al have mutliple hits that can be detected. The total charge for 

1B5 these events however is not significantly different from the events with just 

1B6 one interaction. For this analysis we did not treat them separately since we 

187 were primarily interested in the total charge deposited on the anode. Better 

188 separation of multiple photon interactions on an event by event basis will 

189 be possible with finer segmentation of the readout electrodes and shorter 

190 shaping time. 

191 The purity of the LXe has an impact on charge collection. In the current 

192 setup we achieved an electron lifetime of 200 jis using purification in the gas 

193 phase with heated getter^. We estimated that the level of purity obtained 

194 would result in a loss of 8% of the electrons due to attachment. 

195 For the analysis in this paper we selected events where no net charge 

196 was measured on A2. By demanding the absence of charge on A2 the region 

197 of Al in which events were accepted was smaller than its physical size since 

198 charge depositions close to the edge of Al induced charge on A2. The 

199 effective radius of the tube in which events were accepted was estimated to 

200 be 0.45 cm compared to the Al radius of 0.5 cm. 

201 Figure m shows the distribution of charge due to 511 keV photons in- 

202 cident on the chamber as measured on Al as a function of the drift time 

203 for a 1 kV/cm drift field. The shape of the distribution is the same for all 

204 drift fields. The 511 keV band rises sharply in less than 1 ps, and then falls 

205 slowly until the cutoff which corresponds to the edge of the chamber. The 

206 sharp rise corresponds to photons interacting between the grid and anode. 

207 In that case the electronics, which is not sensitive to the charge induced 

208 by the much slower drifting ions, measures only a fraction of the charge 

209 which is approximately proportional to the distance between the anode and 

210 the interaction point. When the interaction point is between the grid and 

211 the cathode, the measured charge should be independent of the interaction 

212 position. The decline of measured charge with increasing drift times is due 

213 to electron attachment by impurities in the LXe. 

214 Compton scattering interactions are evident below the 511 keV band. 

215 They are due to photons entering the chamber with less than 511 keV 

216 because they have scattered in the passive detector material, mostly the 2 

217 cm of LXe between the vessel wall and the cathode, and to photons escaping 

218 after a Compton scattering interaction in the liquid. 



problem occurred with the purification system during the data taking with the 
Nal coincidence trigger used in this paper resulting in an electron lifetime of 90 jis for 
much of the data presented here. 




Drift time (^is) 

Figure 4: Charge collection as a function of drift time for a 1 kV/cm drift field. The 
curve is a fit based on parametrization obtained from current calculations for energy 
deposits of 511 keV. The scale on the right corresponds to the number of events that 
occurred at a specific time with a certain charge deposition. 



219 We performed a fit of the 511 keV band to extract the drift velocity Vd, 

220 the total charge Qtot produced in the photon interaction and the attenua- 

221 tion length. The values obtained for these quantities are listed in Table [T] 

222 with their statistical uncertainties from the fits. The charge yield Qtot/Qo 

223 is also shown along with Qtot which is the measured charge corrected for 

224 attachment and electronics calibration and Qo is the ratio of the energy 

225 deposited by the 7-ray and the average energy to produce an electron ion 



226 pair: Qo = E^/W with 1^=15.6 eV [17 



227 Figure shows the comparison of our results for the charge yield to 

228 the values obtained in and [iS]. Our results lie in between the two 

229 previous measurements of the charge yield. The obtained drift velocity was 



230 in agreement with previous measurements in [19 



Ed 




Qtot 
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Qtot 


[kV/cm] 


[cm/jis] 


(511 keV e") 


[M 


Qo 


0.33 


0.16 ± 0.01 


19 707 ± 55 


94 ± 3 


0.60 


1 


0.18 ± 0.01 


23 372 ± 59 


61 ± 2 


0.71 


2 


0.20 ± 0.01 


25 092 ± 100 


76 ± 5 


0.77 


2.66 


0.20 ± 0.01 


24 761 ± 35 


60 ± 1 


0.76 



Table 1: Drift velocity (ud), number of electrons (Qtot), electron lifetime r and charge 
yield observed for different electric fields (Ed). 
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231 5. Light Collection 

232 Scintillation light was detected by the LAAPDs located on two sides of 

233 the chamber. Figure shows the sum of the number of photons measured 

234 by both LAAPDs as a function of the electron drift time for the events 

235 where all the charge was collected on the central anode. 

236 The bell shape in Fig. [6] is due to variations of the solid angle with 

237 drift distance which can be calculated by integrating over the LAAPD area 

238 for a given location in the chamber assuming no reflections occurred in 




Figure 6: Light collection as a function of drift time for a 1 kV/cm drift field at 511 keV. 
The curve is a fit based on parametrization obtained from solid angle calculations. The 
scale on the right corresponds to the number of events that occurred at a specific time 
with a certain charge deposition. 



239 the chamber walls. The solid angle varied significantly with the position 

240 of the photon interaction. The arrival time of the electrons provided a 

241 handle on the solid angle variation in the drift direction. However, there 

242 was no information about the position of the interaction within the disk 

243 defined by Al. When the LAAPDs were used independently, the solid 

244 angle variation within this disk introduced a 22% fluctuation in the light 

245 collection. Combining both LAAPDs reduced the fluctuation to 6%. 

24« We fitted the distribution in fig. [6] using a parametrization of the solid 

247 angle, with the total number of photons and the drift velocity as free pa- 

248 rameters. The fit parameters are shown in Table [21 The drift velocity is 

249 consistent with the one extracted from the fit to the charge distribution. 

250 The total number of photons drops with increasing drift voltage in agree- 

251 ment with previous measurements 0, [isj. The number of photons actually 

252 created within the detector was not extracted in this analysis because we 

253 did not measure the photon detection efficiency (PDE) of the LAAPDs. 

254 One of the LAAPDs detected more photons than the other one, which sug- 



255 
256 
257 
258 
259 
260 
261 
262 
263 
264 
265 



gests that there may be a variation of the PDE between LAAPDfl In 
the results given in Table [21 we have assumed 100% PDE for the LAAPD 
that exhibited a more stable operation and scaled the light measured by the 
other LAAPD accordingly introducing a systematic uncertainty because of 
the unknown efficiency (which may be up to 50%). Uncertainties also orig- 
inated from the fact that the ratio between the mean value measured by 
the two LAAPDs varied between data sets by 10%. The light yield was 
computed using a value of 13.8 eV [l7| needed to create one photon at zero 
drift field resulting in 5*0 = 37029 photons. Ntot is the number of measured 
photons corrected for the solid angle of the geometry but not corrected for 
the photo detection efficiency of the LAAPDs. 
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Vd 


cm/ps] 


Ntot (511 keV e-) 


Ntot/S^ 


0.33 


0.15 ± 0.01 


12 161 ± 1269 


0.33 


1 


0.15 


± 0.01 


10 113 ± 1055 


0.27 


2 


0.20 ± 0.01 


9243 ± 964 


0.25 


2.66 


0.21 ± 0.01 


7936 ± 828 


0.21 



Table 2: Electric field {Ed), drift velocity (wd), number of photons {Ntot) and light yield 
observed (see text) for 511 keV photon interactions. 



266 Figu re [7| shows the comparison of our results with values obtained in 

267 [ol and [l8j ^ If the quantum efficiency of the LAAPDs was 60%, which 

268 later results presented here suggest, our results would be in agreement with 

269 previous measurements. 



270 6. Light and Charge Combination 

271 To study the energy resolution we focused on the central region of the 

272 chamber by selecting events with no charge on A2 and choosing a time 

273 window in the drift direction corresponding to 2 mm drift located on the 

274 axis of the LAAPDs where the light collection is maximal as shown in Fig. [6l 

275 The charge signals were corrected for attenuation and the light signals for 

276 the difference between the two LAAPDs and the solid angle dependence in 

277 the drift direction. Resolution results are given as the standard deviation 

278 (a) of a Gaussian distribution unless otherwise stated. Figure E\ shows 
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Figure 7: Light yield relative to the maximum yield measured by different groups at 
different 7-ray energies: this work marked with [6] with □ and [fSj with Q- 



the analysis of a data run at a drift field of 2.66 kV/cm. Evaluating the 
charge and light signals separately gave energy resolutions of 12.1% for 
light and 5.4% for charge by fitting the spectra (shown in the upper plots of 
Fig. [HI) with a sum of two Gaussians and evaluating the mean and width of 
the 511 keV peak. The energy resolution can be improved significantly by 
combining the information from light and charge using the anti-correlation 



2B5 of the two signals [g, ll8|, |20| . The lower left plot of Fig. [8] shows the linear 

2B6 anti-correlation between the light and charge measurement and the axis 

2B7 of the ellipse. Selecting the 511 keV region of the photo-electric-peak the 

288 correlation angle was obtained from a linear fit which provided the axis 

289 of the charge-light ellipse. Projecting the data points along this axis as 

290 described in j6| gave the overall energy resolution. The correlation angle 

291 given here depended on the detector geometry and the efficiency to measure 

292 light and charge separately. The upper left plot shows the charge spectrum 

293 collected on the anode which is equal to a projection of the correlation along 

294 the light axis. In the upper right plot the projection of the correlation 

295 along the charge axis can be seen, giving the spectrum of the collected 



296 light. The lower right plot demonstrates the improved energy resolution 

297 of the combined spectrum when projecting along the correlation axis and 

298 normalizing to the mean charge. The sum of three Gaussians was used as 

299 the fit function to account for the three contributions to the spectrum: the 

300 Compton region (C), the photoelectric peak (P), and scattered events (S) 

301 which lost energy outside the detector, mostly in the LXe before entering 

302 the chamber. 
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Figure 8: The observed charge spectrum (upper left plot), light spectrum (upper right 
plot), correlation between light and charge signals (lower left plot), and combined spec- 
trum using the correlation (lower right plot) for 511 keV photons with a drift field of 
2.66 kV/cm. The data points in the correlation plot (lower left) that are not part of the 
Compton (C) or the photoelectric peak (P) are due to photons that scattered outside the 
detector (S). The linear fit (solid line) giving the axis of the correlation ellipse is depicted 
as well. The fits shown (solid lines) were made with a sum of 3 Gaussians (upper left 
and lower right plot) or 2 Gaussians (upper right plot). 
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Another variable to quantify the anti-correlation between light and charge 



is the correlation coefficient p [22|. Assuming that the probability for a re- 
combining electron-ion pair to produce a scintillation photon is 1 and the 
detector would be able to measure light and charge with 100% efficiency 
and perfect resolution, p should be -1. Deviation from -1 could be due to 
other sources of fluctuations like density fluctuations or delta electrons as 
discussed in 



23 



and the references within. 
Table [3] gives the results of the analysis for different drift fields. The 
best combined energy resolution reached for these data sets was 4.1% at 
2.66 kV/cm drift field (see below). 
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0.33 


13.5 ± 0.2 


7.3 ± 0.5 


4.7 ± 0.1 


56 
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12.2 ± 0.2 


6.0 ± 0.3 


4.3 ± 0.3 


59 


-0.34 


2 


12.8 ± 0.5 


7.0 ± 0.6 


4.8 ± 0.4 


62 


-0.34 


2.66 


12.1 ± 0.1 


5.4 ± 0.2 


4.1 ± 0.1 


58 


-0.26 



Table 3: Energy resolutions (cr) observed at different drift fields for light and charge 
separately and combined result using the correlation. 



313 7. Discussion of Error Sources 

314 In this section we discuss contributions to the energy resolution that 

315 were due to detector inefficiencies or physics constraints like light-charge- 

316 fluctuations. When a photon interacts there is an initially produced number 

317 of ionization charges and scintillation photons which is modified by recom- 

318 bination dependent on the presence of an applied electric field. Table H] 

319 summarizes the variables used in the calculation of error contributions. 





Charge 


Scintillation light 


Initial 
Final 

Measured 


Qm = A Qf 


S^ 

Sf = Si + QiFj.{Ed)Pe^hu 



Table 4: Parameters used in the discussion of energy resolution as described in the text. 

320 Qi,Qf and Qm {Si,Sf and Sm ) are the numbers of initially produced, 

321 post-recombination, and measured charge (light) signals respectively. Fr{Ed) 



322 is the fraction of electron-ion pairs that recombine for a given electric field 

323 Ed, and Pe~^hv = 10] gives the probability for a recombining electron- ion 

324 pair to produce a scintillation photon. Impurities may capture some elec- 

325 trons, which is accounted for by an attenuation parameter A, which depends 

326 on the electron drift distance. The photo-detectors have a photo-detection 

327 efficiency e and cover a fraction of the total solid angle F^. 

328 7. 1 . Charge 

329 The charge resolution is dominated by electronics noise and charge-light 

330 fluctuation. The charge-light fluctuation is expressed by the fluctuation of 

331 the recombination fraction F^, AF^. The charge resolution can then be 

332 written as 

f^QmV _fENCX ^ / AF. y ^ I- A 

\ Qm J \ Qm J \ 1 E^. J Qm 

333 where the first term on the right describes the electronics noise of the am- 

334 plifier, the second term quantifies the light-charge fluctuation, and the third 

335 contribution is the attachment factor which is negligible. AF^ describes the 

336 fluctuation of the recombination and will also occur in the discussion of the 

337 error sources for the light measurement in the next section. Table [5] gives 

338 values for the error contributions to the energy resolution obtained from the 

339 charge measurement. Also shown is the intrinsic energy resolution found by 

340 subtracting the noise of the electronics from the measured resolution. The 

341 only unknown variable contributing to the error of the energy resolution is 

342 AFr which can be calculated once the intrinsic energy resolution is known. 

343 The values for AF,. obtained can also be found in Table [51 
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res. [%] 
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0.33 


7.31 ± 0.54 


5.03 ± 0.04 


5.3 ± 1.5 


3.2 ± 0.9 


1 


6.04 ± 0.33 


4.22 ± 0.03 


4.3 ± 0.9 


3.1 ± 0.7 
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7.00 ± 0.62 


4.29 ± 0.03 


5.5 ± 1.6 


4.2 ± 1.2 


2.66 


5.43 ± 0.17 


3.48 ± 0.03 


4.2 ± 0.4 


3.2 ± 0.3 



Table 5: Contribution of error sources to the energy resolution obtained from the charge 
measurement. 
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1.2. Light 

The fluctuations in the LAAPD specified as the excess noise factor 
F{M) = 2 + kM is dependent on the gain M and affects the resolution. 
For this setup with k = 0.001, F{M) = 2.5. Furthermore, the LAAPD gain 
< 10^ requires that low noise electronics must be used to further amplify 
the signal, which adds electronic noise ENCs- Another source of fluctua- 
tions arises because the solid angle seen by the photo-sensor may vary on 
an event-by-event basis since the solid angle changes with the position of 
the photon interaction within Al. This fluctuation can be corrected for if 
the interaction position is known well from the ionization signal. 

Neglecting other detection fluctuations, the light signal resolution for 
our setup can be written as: 



where the flrst term on the right represents the electronics noise, the 
second term gives the contribution from fluctuations in the LAAPD gain, 
the third term is the fluctuation of the solid angle due to the position of 
the light creation inside the chamber and the fourth term describes the 
light-charge fluctuation. The contribution of the fluctuation in Pe^hv given 
by the last term is negligible or exactly zero if Pe-*hu = 1- Table [6] gives 
values for the error contributions to the measured energy resolution from the 
scintillation light. The solid angle fluctuation amounted to 5.6% and was 
independent of the drift fleld. Also shown is the intrinsic energy resolution 
found when subtracting those error sources due to the detector from the 
measured resolution. AFr was calculated again and can be compared to 
the values obtained from the charge measurement. The values for AFr 
from both the light and charge measurements are in good agreement within 
statistical errors providing a consistency check for the error analysis. 

7.3. Combination 

Combining the light and charge allows improvement to the resolution 
by canceling the fluctuations of Fr by making use of the anti-correlation. 




Ed 


Measured 


Noise 


LAAPD 


Intrinsic 




[kV/cm] 


res. [%] 


[%] 


fluct. [%] 


res. [%] 


[%] 


0.33 


13.5 ± 0.2 


3.3 


0.46 


9.6 ± 0.4 


5.8 ± 1.3 


1 


12.2 ± 0.2 


4.0 


0.55 


6.8 ± 0.8 


3.3 ± 1.0 


2 


12.8 ± 0.5 


5.1 


0.56 


7.1 ± 1.7 


3.1 ± 1.2 


2.6 


12.1 ± 0.1 


4.7 


0.63 


5.5 ± 0.5 


2.5 ± 0.8 



Table 6: Contribution of error sources to the energy resolution obtained using the light 
measurement. 



373 provided the measured charge is corrected for attenuation and the measured 

374 hght for sohd angle and PDE: 



E. = + (3a) 

= Qf + J^ (3b) 



= Q,{l-Fr) + -^ + Q,Fr (3C) 

= Q. + j^ (3d) 

375 

376 where E^ is the energy measured by combining the charge and hght 

377 signals: In eq. [Ha] the light and charge signals were combined; eq. [3b] and 

378 eq. [3c] made use of the formulas in Table [H and in eq. [3d]Fr was eliminated. 

379 The remaining uncertainty in the combined energy resolution is then: 



AE' 



f ENCs 



\ M 



+ F{M)S^ + 



AEnSm 



+ 



+ 



ENCr, 



A 



+ 



FrQ,{l-Pe- 



P. 



(4) 



where the first term on the right hand side (in brackets) originated 
from the light signal contribution and the second term from the charge 
measurement. The last two terms describing the contributions from the 
binomial statistics of A and Pe^hv are negligible. 



384 Table [7] gives the calculated values for the combined energy resolution 

3B5 ^E(./ Ec- The intrinsic combined energy resolution given in the last column 

3B6 of table [7] was obtained by subtracting AEc/Ecin quadrature from the mea- 

3B7 sured values. The solid angle was 12% and the efficiency e for the LAAPDs 

388 was assumed to be 1. 



Ed 


Meas. comb. 


AEJE, 


Intr. comb. 


[kV/cm] 


res [%] 


[%] 


res. [%] 


0.33 


4.74 ± 0.09 


4.0 ± 0.1 


2.5 ± 0.4 


1 


4.31 ± 0.26 


3.7 ± 0.1 


2.3 ± 1.0 


2 


4.78 ± 0.35 


3.8 ± 0.1 


2.9 ± 1.2 


2.66 


4.14 ± 0.10 


3.3 ± 0.1 


2.5 ± 0.4 



Table 7: Contribution of error sources and corrected intrinsic energy resohition for com- 
bined charge- light measurement. 

389 The intrinsic correlation coefficient was calculated but showed a large 

390 uncertainty due to the impact of the uncertainty on AFr. 

391 Another consistency check of the formulas presented here can be made 

392 by extracting the factor ePg^hu in two different ways from the data. We 

393 cannot disentangle the efficiency from the probability but the product can 

394 be obtained from the correlation angle of the 511 keV cloud by writing 

395 C^Yn/-A as a function of Sm/Fn with the slope m describing the axis of the 

396 ellipse: 



= Qf = m-^ + const. (5) 
A Fn 



Ornrr = arctau 



-1 



e—fhu 



The results are listed in Table [HI The mean value over all runs was 
ePe^hu = 0.60 ± 0.03. 

The other method for extracting tPe^hv is to find the slope from the 
plot of the mean values of light vs. charge (coordinates of the center of the 
ellipse) for each electric field setting. This is shown in Fig. [H The linear fit 
to this data gave a slope of ePe-,hu = 0.7 ± 0.3 which agrees with the value 
obtained from the event-by-event method. This result would be consistent 
with Pe^hv = 1 if the efficiency of the APDs were 60%. 



Ed [kV/cm 




^corr 




0.33 


56.2 


0.67 ± 0.03 


1 


58.7 


0.61 ± 0.11 


2 


62.3 


0.53 ± 0.30 


2.66 


58.5 


0.61 ± 0.05 



Table 8: Calculating ePe^hiy on event by event basis. 
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Figure 9: Method of extracting ePe^hv from mean values of the correlation cloud pa- 
rameters described in the text. The solid line is the linear fit to the data. 




7.4- Improving Light Resolution using Position Reconstruction 

In the current system, the energy resolution contribution from light was 
partially limited by the uncertainty in the position of the light source due to 
the 1 cm dia. size of Al i.e. from the fluctuation of the solid angle within Al. 
Applying the formulas presented above to a Geant4 2J] simulation, it was 
found that knowing the position of the light signal to 1 mm would improve 
the light resolution by about 1.5% and the combined resolution by up to 
0.5% giving 10% for light alone and 3.6% for the combination of charge and 
light, consistent with the resolution reported in iej. Tabled summarizes the 



414 values obtained with the simulation comparing the two cases of not knowing 

415 the position of the interaction within Al and being able to correct for it. 





Light res. [%] 


Charge res. [%] 


Combined res. [%] 


Measured 
Simulated 
Corrected 


12.1 ± 0.1 
12.0 ± 0.2 
10.4 ± 0.2 


5.4 ± 0.2 
5.4 ± 0.1 
5.4 ± 0.1 


4.1 ± 0.1 

3.9 ± 0.1 
3.6 ± 0.1 



Table 9: Energy resolutions obtained from the simulation with (corrected) and without 
(simulated) the correction for the position of the interaction within Al in comparison 
with the measured resolutions. The charge resolution was not affected by the correction. 



416 8. Summary and Conclusion 

417 Measurements have been made of the response of a liquid xenon drift 

418 chamber to irradiation by 511 keV photons. Using a model accounting 

419 for the sources of uncertainty in the energy resolution we also determined 

420 values for the intrinsic energy resolutions. Figure fTOl summarizes the results 

421 for charge, light and combined energy resolution as a function of the drift 

422 field. The main error contribution to the combined energy resolution, apart 

423 from the solid angle fluctuations which can be eliminated by utilizing the 

424 position measurement, originated from the APD gain fluctuation and the 

425 anode noise. Both were about of 2.7%. 

426 Figure fTTl depicts the values for the intrinsic energy resolutions obtained 

427 by subtracting the detector contributions from the values in fig. [TOl The 

428 error bars given are statistical. 

429 Based on these results, the combined energy resolution of <3.5% (or 

430 < 8% FWHM) would be anticipated in a detector configuration suitable for 

431 applications to PET which would have comparable light collection efficiency 

432 to the prototype detector described above and <1 mm spatial resolution. 

433 Reducing the anode to grid spacing to 1 mm and the grid wire spacing to 

434 1 mm will reduce the width of the pulses and minimize the dependence of 

435 the pulse shape on the location of the electron cloud. Further improvements 

436 are foreseen in areas including purification and low noise electronics. 
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Figure 10: Energy resolution from charge (Q) ^^^d light (□) measurements as well as the 
combined (★) resolution for different drift fields. 
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Figure 11: Intrinsic energy resolution from charge (Q) said light (□) measurements as 
well as the combined (★) intrinsic resolution for different drift fields. 
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